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Abstract 

> 

■ We present results from the XMM-Newton observation of the non-coohng flow 

cluster A 1060. Large effective area of XMM-Newton enables us to investigate the 
O ! nature of this cluster in unprecedented detail. From the observed surface brightness 

^ I distribution, we have found that the gravitational mass distribution is well described 

by the NFW profile but with a central density slope of ~1.5. We have undoubt- 
edly detected a radial temperature decrease of as large as ~30% from the center to 
the outer region (r ~ 13'), which seems much larger than that expected from the 
^ ■ temperature profile averaged over nearby clusters. We have established that the tem- 

perature of the region ~ 7' southeast of the center is higher than the azimuthally 
averaged temperature of the same radius by ~20%. Since the pressure of this region 
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^ i already reaches equilibrium with the environment, the temperature structure can be 

interpreted as having been produced between 4 x 10^ yr (the sound-crossing time) and 
3 X 10^ yr (the thermal conduction time) ago. We have found that the high-metallicity 
blob located at ~ 1.'5 northeast of NGC 3311 is more extended and its iron mass of 
1.9 X 10^ Mq is larger by an order of magnitude than estimated from our Chandra 
observation. The amount of iron can still be considered as being injected solely from 
the elliptical galaxy NGC 3311. 

Key words: galaxies: clusters: individual (Abell 1060) - galaxies: intergalactic 
medium - X-ray: galaxies 
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1. Introduction 

Clusters of galaxies, being the largest virialized systems in the universe, are filled with 
the intracluster medium(ICM), which consists of X-ray emitting hot plasma with a typical 
temperature of a few xlO^ K. The ICM gives us rich information about the structure and 
evolution of clusters of galaxies. X-ray Spectra of ICM immediately reveal its temperature 
and metal abundance, and the ICM surface brightness shows us ICM density distribution. 
Through the assumption of hydrostatic equilibrium, furthermore, it is possible to measure the 
distribution of the gravitational mass in clusters. 

The cluster mass distribution is an important observational clue in determining the 
nature of the dark matter and in constraining the structure formation scenario. In particular, 
the CDM model predicts a significant central cusp in the dark matter profile 

(Navarro, Frenk, & White 1996), in contrast to the fiat core profile described by the /3 
model. High- sensitivity X-ray study is a powerful method to look into the gravitational mass 
distribution (e.g. Churazov et al. 2003). For clusters characterized by cool central cores, pres- 
ence of the bright central galaxy greatly hampers the mass determination for the pure cluster 
component. Therefore, detailed study of a non cool-core cluster has a significant importance. 

Recent observational studies of the clusters with Chandra and XMM-Newton, with their 
powerful imaging capability and large effective area, unveiled a lot of new interesting features. 
Mazzotta et al. (2003) found an unusual X-ray morphology, blob and hole like features, in 
the central dense region of the 2A 0335-1-096 cluster. Markevitch et al. (2000) detected sharp 
brightness edges, called cold fronts, with lower temperature inside the subcluster in A 2142 
cluster. Furusho et al. (2003) found two high-metallicity blobs located symmetrically with 
respect to the center of the poor cluster AWM 7. Piffaretti et al. (2005) investigated temperature 
profiles of 13 nearby cooling fiow clusters observed with XMM-Newton, and found that the 
temperature decreases by ~30% between 0.1 r^j^ and 0.5 Vmr- This feature is in good agreement 
with the result of Markevitch et al. (1998) who found that the temperature of 30 nearby clusters 
show a significant decline at large radii. Recent Chandra observation of A 1060 revealed that 
the central two elliptical galaxies emit X-ray (Yamasaki et al. 2002), and the temperature and 
the abundance distribution are unexpectedly inhomogeneous (Hayakawa et al. 2004). 

A 1060 is a nearby cluster [z = 0.0114) characterized by a smooth and symmetric ICM 
distribution. The ASCA and ROSAT observations found that the temperature and the abun- 
dance showed very uniform distribution (Tamura et al. 1996; Furusho et al. 2001). For these 
reasons, this cluster was thought to be in a highly relaxed state, and has been used to estimate 
the gravitational mass profile which can be directly compared with the results of hydrodynamic 
cluster simulation (Tamura et al. 2000). 

In this paper, we report the results of XMM-Newton observation of the A 1060 cluster. 
Based on the surface brightness distribution, we have computed the total gravitational mass of 
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the cluster. We confirmed the existence of the high-met allicity blob discovered by our Chandra 
observation but with a larger iron mass content. Using the high spatial resolution and the 
large effective area, we searched for a sign of a sub-cluster merger from a two-dimensional 
temperature map. Throughout this paper, we assume Hq = 75km s~^ Mpc~^ with go = 0.5. 
Accordingly, an angular size of 1' corresponds to 13 kpc. The solar number abundance of iron 
relative to H is assumed to be 4.68 x 10~^ (Anders & Grevesse 1989). We employ Galactic 
absorption as Nu = 4.9 x 10^°cm~^. 

2. Observation and Data reduction 

The XMM-Newton observation of the A 1060 cluster of galaxies was carried out on 2004 
June 29 for a total exposure time of 64 ksec. All the EPIC cameras were operated in the 
Full frame mode with the medium filter inserted. Data reduction and analysis were performed 
with SAS version 6.0 and ClAO 3.0. Unfortunately, part of the data was affected by high 
background flares. In order to remove them, we adopted the following data cleaning process. 
We made light curves in the bands 10-12 keV for MOS and 12-14 keV for pn, where the signals 
are dominated by background events, and then we removed the time interval when the count 
rate is larger than 0.3 c s"^ for MOS and 0.8 c s~^ for pn. Then, we selected the photon events 
with pattern 0-12 for MOS and 0-4 for pn, and flag=0 ^ . After these selections, some 35.9, 
36.4, and 25.2 ksec of good time intervals remained for MOSl, M0S2, and pn, respectively. 

The blank sky data sets of Read & Ponman (2003) were used to subtract the background 
of our data sets. We extracted the background events using the same selection criteria as above. 
The source-to-background count-rate ratio was calculated from the count rates in the bands 
10-12 and 12-14 keV for MOS and pn (see Katayama et al. 2004), which resuh in 1.132, 1.250, 
and 1.200 for MOSl, M0S2, and pn, respectively. After subtracting the blank sky data, being 
scaled by these factors, from the source data sets, we created X-ray images for each detector 
separately. Figure 1 (a) and (b) show the X-ray image of MOSl in the band 0.8-8.0 keV. The 
telescope vignetting effect had been corrected using the exposure map which was calculated 
with the SAS "eexpmap" ^ command. The image had been smoothed by a Gaussian with a la 
width of 4.4" using the ClAO csmooth command. The two bright elliptical galaxies (NGC 3311 
and NGC 3309) are clearly visible in the central region of the cluster (Fig. 1 (b)). Properties 
of these galaxies have been reported by Yamasaki et al. (2002). A blob-like high brightness 
region was located at northeast of NGC 3311. Except for these central features, the whole 
surface brightness distribution remains nearly spherically symmetric. Although the surface 
brightness map in Fig. 1(a) seems elongated in the north-south direction, the lengths of the 
semi-major /minor axes of the isophote at r ~ 5', for example, agree within a few percents. 

^ http:/ /xmm. vilspa.esa.es/sas/current/doc/evselect/ 
^ littp:/ /xmm. vilspa.esa.es/sas/current/doc/eexpmap/ 
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Fig. 1. (a) Adaptively smoothed MOSl X-ray image of the 30' x 30' region, overlaid with the inten- 
sity contours. The contours were logarithmically scaled and divided into 20 spaces between 10^^ and 
10"'^ counts s~^ range. The red circle represents XMM-Newton field of view, (b) Close up view in the 
central 8' x 8' region. 

We created on-source spectra of each detector within a radius of 13' from the cluster cen- 
ter. The point-Hke sources seen in Fig. 1 were removed using SAS "edetect_chain" ^ command. 
The central two elliptical galaxies possess an extended halo (Yamasaki et al. 2002). Accordingly, 
we have excluded manually the circular regions with radii of 25" and 20" centered on NGC 3311 
and NGC 3309, respectively, for accumulating spectra in order to avoid contamination from the 
halo. The blank sky spectra were subtracted from the on-source spectra after the appropriate 
scalings to remove the background photons. For pn data, there are strong fluorescence lines of 
Ni, Cu and Zn at around 8 keV. Since these lines can not be completely removed by the back- 
ground subtraction, we introduce Gaussians in the spectral model to represent these lines. First 
of all, we fitted the spectra from each detector separately with the absorbed MEKAL model 
in the bands 0.8-8.0 keV for MOSl, 2 and 0.8-10.0 keV for pn. As summarized in table 1, with 
the absorption fixed at the galactic value of A'^h = 4.9 x 10^° cm"^ (Dickey & Lockman 1990), 
the temperature is distributed in the range 3.10-3.30 keV among the three detectors. If the A^h 
value is set free to vary, a similar temperature distribution {kT = 3.13-3.31 keV) is obtained. 
The Nu value from the pn fit, however, shifted to a lower value of 1.2to;4 x 10^° cm~^, which is 
significantly smaller than 4.9 x lO^'' cm~^ obtained from the radio 21 cm measurement (Dickey 
& Lockman 1990). We then carried out a combined fit of the absorbed MEKAL model to 
the MOSl-l-2 and pn spectra of the r < 13' region in the band 2.0-10.0 keV, free from the 

^ http:/ /xmm. vilspa.esa.es/sas/current/doc/edetectjchain/ 
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Table 1. Best-fit Parameters of the MEKAL model to the Spectra in the central r < 13' region 



Inst. 


kT 


Z 




xVd.o.f 




[keV] 


[Zq] 


[xl020 cm-2] 




MOSl 


o on+0.04 
'-'■'-"-'-0.04 


41+0 02 

^■^^-0.02 


4.9(fix) 


1.244 




o 9o-f0.05 
''■^'-'-0.05 


"■^"-0.02 


"•"-0.6 


1.226 


M0S2 


o 90+0.03 
'^■^"-0.04 


41+0 02 

"■^^-0.02 


4.9(fix) 


1.366 




q 1 q + 0.05 

''■^'^-0.05 


40+0 02 

^■^'-'-0.02 


6 1+0-^ 

"■^-0.6 


1.344 


pn 


o 1Q+0.03 
'-'•^'-'-0.02 


^■'J'-o.oi 


4.9(fix) 


1.475 




n qi +0.02 

'-'•'->^-0.02 


sq+0 02 

^■'-'^-0.02 


-1 9+0.5 
-■^•^-0.4 


1.238 


MOS+pn 


> 2 keV 


o oq+0.04 
'^■'-'^-0.04 


n S'S+0 02 
""^^-0.02 


4.9(fix) 


1.145 


> 0.8 kcY^ 


q 97+0.02 
"^■^ ' -0.02 


'?q+o.02 

^■'-'''-0.01 


4.9(fix) 


1.250 



Note — All errors are at the 90% confidence level. 

a: The energy bands for MOS and pn are 0.8-8 keV and 1.6-10 keV, 

respectively. 
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Fig. 2. MOSl+2 and pn spectra extracted fi'om the r< 13' region fitted in the bands 0.5-8 keV and 
0.5-10 keV, respectively, with the MEKAL model undergoing the galactic absorption. 

absorption, and simply retrieved the data points below 2 keV. The result is shown in Fig. 2. 
The hydrogen column density is fixed at the galactic value, 4.9 x 10^'' cm~^. The best-fit pa- 
rameters are summarized in table 1. The MOS and pn results are inconsistent below ~1 keV, 
and the pn fit requires a strong soft excess. A similar phenomenon has been reported from the 
XMM-Newton observation of A 1413 (Pratt & Arnaud 2002). According to these authors, this 
effect could be due to a true soft excess emission (e.g. Durret et al. 2002) and/or an artifact 
due to remaining calibration uncertainties. In the case of A1060, however, such a soft excess 
is incompatible with previous observations (Tamura et al. 2000). Accordingly, we neglect the 
energy band below 1.6 keV for the pn, and hereafter adopt the 0.8-8.0 keV and 1.6-10.0 keV 

5 



I 

i 



•3 



•1 

o 
•d 



0.6 



0.8 - 



lOOkpc 



SOOl 



-t — I — I — I — I — I — I — I- 



O MOS{0.e-B.Ok#V)+pA 
-I- Cliandra(O.9-B.0keV] 



H — \ — I — I — I — I — I — } — I — 1 — h 



— I I I ■ I I I I i_ 



10 15 
Ra<Uus[min] 



20 



Fig. 3. Radial distributions of tlic temperature (top), metal abundance (middle), and reduced value 
(bottom) of the spectral fits to the concentric annular regions. The red diamonds show the results from 
the MOS spectra in the bands 0.8-8.0 keV (pn energy band is fixed in 0.8-8.0 keV). The black dots show 
the result of Chandra observation. The error bars are at the single-parameter 90% confidence level. 

bands for the MOS and the pn, respectively. The results of the combined fit in these energy 
bands are also summarized in table 1, showing that the temperature of the plasma within the 



-1-0.02 y 
-0.01^0- 



projected central r < 13' region is 3.27±0.02 keV with the metal abundance of 0.39 
3. Temperature and Abundance Distribution 



3.1. Radial Profile 

We have produced radial temperature and abundance profiles in annuli centered on the 
position of NGC 3311. The central two elliptical galaxies (NGC 3311, NGC 3309) and the 
other point-like sources are removed (§ 2). The annuli are segmented from inner to outer radii 
so that they include at least 20,000 photons for all the detectors. The resulting widths of the 
annuli are ~30" or larger, which is larger than the HPD of the XRT (~17"). The MOS and pn 
spectra in the bands 0.8-8.0 keV and 1.6-10.0 keV, respectively, are fitted simultaneously with 
an absorbed MEKAL model. The absorber's hydrogen column density is fixed at the Galactic 
value of 4.9 x 10^° cm~^. Figure 3 shows the result of radial temperature and abundance 
distributions. For comparison, we have overlaid the Chandra results (Hayakawa et al. 2004). 
Owing to large effective area of XMM, the gradual temperature decline from kT = 3.5 keV to 
2.5 keV toward the outer radii, which is not very clear in the Chandra result, is unambiguously 
detected. Because the field of view of the Chandra observation was offset toward the southwest 
direction, the outer two annuli were of poor statistics and probably affected by background 
uncertainty. The radial abundance increase of ~20 % in the inner 3' region also becomes very 
clear for the first time. 

We have so far adopted the position of NGC 3311 as the cluster center simply because 
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it locates at the apparent peak of the surface brightness, and one may suspect that the radial 
profiles shown in Fig. 3 could be significantly different if we draw them centered on a different 
position. Accordingly, as an alternative, we have tried to define a cluster center with the 
geometry of the isophotes. We have fitted a circle to every isophote shown in Fig. 1(a) with its 
central position and radius being set free to vary, and have found that the center of the circle 
fitted to the second largest isophote (r ~ 11') shows the largest shift of 22" from NGC 3311 
in the western direction. We have made the radial profiles again with this new center. The 
resultant difference from the previous results is, however, quite small, which is only 5% and 
15% at most for the temperature and the abundance over the entire radius range, respectively. 
In particular, the 5% random difference in the temperature is much smaller than the systematic 
temperature drop from 3.5 keV at the center to 2.5 keV at r = 13'. Accordingly, we adopt the 
position of the galaxy NGC 3311 as the cluster center throughout this paper. 

3.2. Two-Dimensional Distribution 

3.2.1. Temperature distribution 

By taking full advantage of large effective area of XMM-Newton, we have attempted to 
draw a two-dimensional temperature map on the basis of spectral analysis. We divided each 
annulus which was adopted to obtain the radial profiles in Fig. 3 evenly into 8 azimuthal sectors. 
The radial width of annuli was determined by requiring a minimum of 8000 counts for all the 
eight sectors. The spectra were fitted, after background subtraction, with a single-temperature 
MEKAL model with the galactic absorption. This process is the same as that described in 3.1. 
Figure 4(a) and (b) show the resulting two-dimensional temperature distribution, and its radial 
profiles in the eight azimuthal directions, respectively. For reference, the azimuthally-averaged 
radial temperature profile (Fig.3) is overlaid in panel (b). It is remarkable that the r = 4' — 8' 
region in the two south-eastern sectors (labelled SE and ES in the figure) seems to be hotter 
than the other sectors at the same radii. 

3.2.2. Metal distribution 

We reported several high-met allicity regions near the center of A 1060 cluster (r < 5') 
based on our Chandra observation (Hayakawa et al. 2004). We thus attempted to produce fully 
two-dimensional metal abundance map in the central r ^ 4' region. We divided the sky region 
into rectangular cells with 22" on a side, and evaluated the metal abundance by fitting MEKAL 
model to the spectra obtained from 3x3 cells. We have obtained two-dimensional metallicity 
distribution through taking running average by shifting the integration region stepped by one 
cell. In order to keep statistical quality, we discard any 3x3 cell if it does not contain at least 
4000 counts. The result is shown in Figure 5(a). A typical statistical error of the metallicity is 
10-15%. We note that there is clearly a high-metallicity region at ~1.'5 northeast of NGC 3311 
whose position is consistent with the blob detected with the C/iandra observation. The projected 
metallicity within a radius of 1' centered on the blob is 0.65 ±0.06^0, which is compatible with 
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Fig. 4. (a): Temperature map obtained from speetral fits. The overlaid contours show smoothed X-ray sur- 
face brightness. The three hues indicate MOSl(red), M0S2 (green), and pn(blue) fields of view. (b):Radial 
temperature profiles of the 8 sectors shown in (a). The red line shows the azimuthally-averaged radial 
temperature profile shown in Fig. 3 (top panel), ignoring errors. 

the Chandra value of O.TAIIq^^Zq, and is about twice as large as in other regions. Figure 
5(b) shows a radial metallicity distribution in the concentric annuli centered on the blob. 
Assuming that the metallicity of the ICM out of the blob region is uniform at ~ OAZq, the 
deprojected central metallicity of the blob is obtained to be ^1.5Zq, which is consistent with 
the Chandra value (Hayakawa et al. 2004). Note, however, that the radial extent of the high- 
metallicity region is found as large as r ~ 2' (Fig. 5b), whereas it was only ~40" in the Chandra 
observation. This is because the brightness of the blob is low and close to the photon limit 
of the Chandra ACIS. None of the other high- metallicity blob-like regions suggested by the 
Chandra observation, which are all fainter than the central blob, was significantly detected in 
the XMM-Newton observation. 



4. Discussion 



The present XMM-Newton observation of A 1060 provided several pieces of new knowl- 
edge on the surface brightness, temperature and abundance distributions in unprecedented 
detail. They include the significant drop of the temperature toward outer radii (§ 3.1), the 
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Fig. 5. (a) Metal abundance distribution in the central r < 5' region based on spectral analysis. The pixel 
size is 22" x 22" (5kpc across). Spectral analysis is carried out for a square region of 3 x 3 pixels. The 
overlaid contours show smoothed X-ray surface brightness, (b) Radial metallicity distribution in annuli 
centered on the position of the high- metallicity blob, according to the concentric annuli of (a). 

existence of the south-eastern high temperature region (§ 3.2.1), the high-metalhcity blob in 
the central region (§ 3.2.2), and so on. We discuss their implication in some detail. 

4.1. Mass Distribution 

We first derive the gravitational mass profile based on the density distribution according 
to the following steps (see also § 4.2 of Hayakawa et al. 2004): 

1. We first calculate the radial surface brightness profile from the MOSl-l-2 image in the 
0.8-8.0 keV band, after application of the SAS task "evigweight" in order to eliminate the 
telescope vignetting effect. 

2. Under the assumption of constant temperature and spherically symmetric distribution 
of the ICM, we fit the radial surface brightness profile in the range r > 3 with a single 
/9-model. The fitted parameters are = 7'.3 ± 1^8 and /3 = 0.691^;^^, respectively. 

3. The best-fit /5-model is extrapolated to the center, and the ratio of the data to the model 
is calculated as a function of the radius. The radial profile of the surface brightness ratio 
is fitted with an exponential function. 

4. The density of the plasma p(r) predicted by the /5-model is enhanced by multiplying the 
square root of the above ratio, and the new radial profile of the brightness is calculated. 

5. The procedures 3 and 4 are repeated until the residual between the data and the model 
practically vanishes. 
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Fig. 6. (a) Gravitational mass profile estimated from the single /3- model (upper dashed line) and the one 
obtained from our analysis (upper solid line) under the assumption of the isothermal ICM in hydrostatic 
equilibrium. The lower lines represent ICM masses only. The radial temperature gradient (Fig. 3) is 
reflected for the diamonds, (b) Gravitational mass density profile based on (a). The dashed-dotted line 
shows the Moore's model. Two representative slopes, r~^-^ and r^^, are also shown by solid lines. 

As a result, we obtained the central electron number density to be uq = 11.7~^q'1 x 10"^ cm~^. 
This value is slightly larger than that obtained from the Chandra observation uq = 8.2jli;o ^ 
10"'^ cm~^ (Hayakawa et al. 2004). This is mainly because we have substituted the exponential 
function for the Gaussian adopted for the Chandra data in the procedure 3 above. Since the 
statistics of the radial surface brightness profile of Chandra data was poor at the cluster center, 
the choice of the model to fit to the data-to-model ratio, especially at the cluster center, is 
not unique. We thus simply adopted the Gaussian with no particular reason. In analysing 
the high quality XMM-Newton data, however, we have found that the exponential function is 
significantly better than the Gaussian. 

Given the radial density profile that can reproduce the surface brightness distribution, 
we then calculate the total gravitational mass with the equation 

rkT{r) f d\np{r) d\nT{r)\ 



tot 



HrripG \ dliar d\nr J 
under to assumption of the hydrostatic equilibrium and the spherical symmetry of ICM. The 
integrated gravitational and gas mass profiles are shown by the solid lines in Fig. 6(a) under 
the assumption of the constant temperature T(r) = (T) (=3.3 keV; § 2). The corresponding 
mass density profiles are also shown in Fig. 6(b). For comparison, those assuming the single 
/9-model are plotted with the red dashed lines. Shown by diamonds are the profiles that reflect 
the radial temperature variation (Fig. 3). Owing to large effective area of XMM-Newton, we 
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have unveiled the significant concentration of mass in the central region that can be represented 
by a density cusp with oc r~^'^. In order to confirm this, we tried a brief model fitting to the 
observed profile in Fig. 6(b). We tried two analytical models; the NFW model (Navarro, Frenk, 
& White 1997) and the Moore's model (Moore et al. 1999). They have a central density cusp 
proportional to and r~^'^, respectively. As a result, we found that Moore's profile could be 
better fit to the data, which is displayed in Fig. 6(b) with the dash-dotted line. The central cusp 
of the gravitational mass density profile in A 1060 was first suggested by Tamura et al. (2000) 
based on the analysis of the ASCA and ROSAT data.. Their analysis, however, did not take into 
account the compact halo associated with the elliptical galaxy NGC 3311 which could raise the 
cusp-like density enhancement. In the analysis of Chandra data, we could remove contribution 
from NGC 3311. Poor statistics, however, precluded us from detecting the density cusp of the 
total mass (Hayakawa et al. 2004). 

The gravitational mass profile allows us to compute ta for an overdensity A using 

where pc is the critical density (Piffaretti et al. 2005). In the SCDM model, the virial radius 
corresponds to A = 180. Our total mass distribution indicates rigo to be about 1.35 Mpc. This 
is in good agreement with the theoretical prediction by Evrard, Metzler, & Navarro (1996), 
ri8o = 2.6 Mpc {Tx/lO keV)^/^ giving rigo = 1-49 Mpc for Tx = (T) = 3.3 keV, as the mean 
emission-weighted temperature. We hereafter adopt rigo = 1.35 Mpc. 

4-2. Abundance Distribution 

The metal abundance distribution in the central r < 5' region (Fig. 5) indicates a high- 
metallicity region centered at ~l!5 northeast of NGC 3311. This enhancement was identified 
in our Chandra paper as the metal rich blob with a radius of 40" (= 8.7 kpc). We estimated 
the total mass of iron contained in the blob to be 3.0 x 10^ Mq. XMM-Newton data showed, 
however, that the size of the high-met allicity region is significantly more extended than the 
Chandra result even considering an image blur of the XMM-Newton data. This is probably 
due to the photon limit of Chandra data at the iron line energy band 6-7 keV. We thus re- 
estimated the mass of iron contained within a radius of 1.'5 (~ 20 kpc) from the blob center. As 
we already mentioned in § 3.2.2, the projected iron abundance measured with XMM-Newton is 
0.65±0.06Zq, which is consistent with the Chandra value of O.TA^q'^IZq. After the deprojection 
as we did for the Chandra data (Hayakawa et al. 2004), the total mass of iron in the blob is 
estimated to be 1.9 x 10^ Mq, which now becomes larger roughly by an order of magnitude. 

The stellar mass of NGC 3311 is estimated to be 4.8 x 10^^ Mq from the relation M/Lb = 
8.5{Mq/ Lq) with niB = 11.15 (Faber et al. 1989). Assuming that the average iron abundance 
in the stellar mass of NGC 3311 is IZq, we obtain the total stellar iron mass as ~1.3 xlO^ Mq. 
As a result, the iron mass of the northeast blob is ~ 1.5% of the stellar iron content, and 
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can reasonably be supplied by the single galaxy NGC 3311. Schindler et al. (2005) simulated 
metallicity evolution in galaxy clusters which includes interaction between the galaxies and 
the ICM. They have noted that the ram-pressure stripping is more efficient than the galactic 
wind in transferring metal-enriched interstellar medium into the ICM, and the expelled metals 
are not mixed immediately with the ICM. The northeast high-metallicity region, likely to be 
associated with NGC 3311, is reminiscent of the ram-pressure stripping from the X-ray image. 
If this is the case, NGC 3311 is moving with a velocity larger than a few 100 km~^ s in the 
central region of A 1060. The ram-pressure stripping is also likely to be the reason why this 
central galaxy shows a very compact X-ray halo (Yamasaki et al. 2002). 

4-3. Temperature Distribution 

Plasma temperature of the A 1060 cluster decreases from the center to r = 13' by ~30 %, 
which is a larger decrement than that detected in our Chandra observation (~20 %). This 
is probably because background uncertainty of the Chandra observation in the outer r ^ 10' 
region. Our result is consistent with that of Tamura et al. (2000) claiming that the temperature 
decreases down to about 2.4 keV at a radius of ~ 12' — 20' based on their analysis of the ASCA 
and ROSAT data.. 

Markevitch et al. (1998) observed 30 nearby clusters and obtained remarkably similar 
trend of the temperature normalized by the emission weighted average temperature (T), as a 
function of the radius normalized by rigQ. We thus have attempted to plot our temperature 
profile with rigo = 1.35 Mpc (§ 4.1) and (T) = 3.3 keV (§ 2). The resulting temperature profile 
is shown by diamonds in Fig. 7, while the composite profile of Markevitch et al. (1998) is shown 
as the region surrounded by the dashed line. The profile of A 1060 obviously becomes out of 
the allowed range in the region r/rigo ^ 0.07. One possible reason to cause such a deviation 
may be an extra heating in the central region, due to a sub-cluster merger, for example. If we 
are allowed to take (T) = 2.5 keV, for example, the normalized temperature profile well fits in 
the Markevitch's range, as displayed in Fig. 7. 

We also recognize that there is a high temperature region at r ~ 7' from NGC 3311 in the 
southeast direction (Fig. 4). This feature can be confirmed also by the color-coded temperature 
map based on the HR analysis in Furusho et al. (2001). Figure 8 shows the radial surface 
brightness profile of the southeast direction, compared with that in the northwest direction in 
the top panel. The ratio of the two profiles (middle panel), and the radial temperature profile 
in the southeast direction (bottom panel) are also shown. The surface brightness around the 
position of the southeast hotter region is slightly fainter than that in the northwest region. The 
pressure P and the emissivity I are represented by the electron density rig and the temperature T 
as P oc UeT and / oc ngT^/^. Combining these relations by eliminating rie, we obtain P oc p/'^T^/^. 
At the southeast high temperature region, the surface brightness is lower than the northwest by 
~ 25%, while the temperature exceeds the hnearly-fitted value by ~ 20%. Thus, the southeast 



12 




Fig. 7. Temperature distribution as a function of the radius from the cluster center in a unit of rigo- The 
diamonds and sohd hue correspond to the red diamonds in Fig. 3, normahzed by the emission- weighted 
average temperature < T > of 3.3 keV and 2.5 keV, respectively. The dashed line define a region that 
encloses most temperature profile and their error bars of Markevitch et al. (1998). 
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Fig. 8. Radial surface brightness profile in the SouthEast direction in comparison with the NorthWest 
(NWired) (top), the ratio of these two profiles (middle), and the temperature distribution in SE direction 
(bottom). The dashed line in the bottom panel is the linear fit to the data in the range r < 3' and r > 10'. 
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region is in pressure equilibrium with the environment. The diffusion time of pressure gradient 
is of the same order as the sound-crossing time. The time for reaching pressure equihbrium 
over the scale 3' (~ 70 kpc) can be estimated to be about 4 x 10^ yr, for the sound speed of 
~ 1000 km s~^. On the other hand, the temperature difference of 1 keV over 3' distance is 
relaxed by thermal conduction whose time scale is estimated to be about 3 x 10® yr in the case 
of A 1060. It is likely that the temperature structure has been produced between 3 x 10® yr 
and 4 x 10^ yr in the past. 

5. Conclusion 

We have presented the XMM-Newton observation of the A 1060 cluster of galaxies. 
Large effective area of XMM-Newton enables us to investigate spatial distribution of the surface 
brightness, the temperature and the metal abundance in unprecedented detail. We have derived 
a gravitational mass profile based on a model-independent estimation of the density profile. The 
overall shape can be well represented by Moore's profile with a central density slope oc r~^'^. 
The temperature is found to decrease undoubtedly by 30 % from the center to the edge of the 
field of view r = 13'. This drop is clearly steeper than the Markevitch's composite profile. 

Our spatially-resolved analysis revealed slight temperature enhancement at r ~ 7' from 
NGC 3311 in the southeast direction. It is thought that some heating, like a subcluster merger, 
has occurred in this region in the past between the sound crossing time 4 x 10^ yr and the 
thermal conduction time 3 x 10® yr ago. We note that the deviation of temperature profile of 
A 1060 from the Markevitch's one can also be explained if the central region is heated like the 
way seen in the southeast region. 

The high resolution metallicity map for the central r < 5' region confirmed the high- 
metallicity blob discovered with the Chandra observation locating at around 1'5 northeast of 
NGC 3311. Its spatial extent is, however, found to be as large as r ~ 2', which is significantly 
larger than that estimated with Chandra (r ~ 40"). This discrepancy is probably because the 
Chandra observation is photon-limited in the iron line band. The iron mass in this region is 
revised to be 1.9 x 10^ Mq. Although this is larger than the Chandra value roughly by an order 
of magnitude, it can still be supplied solely by NGC 3311. 
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